Levels of Complexity in Scale-Invariant Neural Signals 
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Many physical and physiological signals exhibit complex scale-invariant features characterized by 1// scal- 
ing and long-range power-law correlations, suggesting a possibly common control mechanism. Specifically, it 
has been suggested that dynamical processes influenced by inputs and feedback on multiple time scales may 
be sufficient to give rise to 1/f scaling and scale invariance. Two examples of physiologic signals that are 
the output of hierarchical, multi-scale physiologic systems under neural control are the human heartbeat and 
human gait. Here we show that while both cardiac interbeat interval and gait interstride interval time series 
under healthy conditions have comparable 1/f scaling, they still may belong to different complexity classes. 
Our analysis of the magnitude series correlations and multifractal scaling exponents of the fluctuations in these 
two signals demonstrates that in contrast with the nonlinear multifractal behavior found in healthy heartbeat 
dynamics, gait time series exhibit less complex, close to monofractal behavior and a low degree of nonlinearity. 
These findings underscore the limitations of traditional two-point correlation methods in fully characterizing 
physiologic and physical dynamics. In addition, these results suggest that different mechanisms of control may 
be responsible for varying levels of complexity observed in physiological systems under neural regulation and 
in physical systems that possess similar 1/f scaling. 

PACS numbers: 05.40+j, 05.45Tr, 87.10.+e, 87.19.Hh, 87.45Dr., 87.23Ge, 87.80.-y, 87.90. +y 



I. INTRODUCTION 

Many dynamic systems generate outputs with fluctuations 
characterized by 1/ /-like scaling of the power spectra, S(f), 
where / is the frequency. These fluctuations are often asso- 
ciated with nonequilibrium dynamic systems possessing mul- 
tiple degrees of freedom [1, 2], rather than being the output 
of a classic "homeostatic" process [3-5]. It is generally as- 
sumed that the presence of many components interacting over 
a wide range of time or space scales could be the reason for the 
1/f spectrum in the fluctuations [6, 7]. Fluctuations exhibit- 
ing l//-like behavior are often termed "complex", since they 
obey a scaling law indicating a hierarchical fractal organiza- 
tion of their frequency (time scale) components rather than 
being dominated by a single frequency. 1/f behavior is com- 
mon in a variety of physical, biological and social systems [7- 
15]. The ubiquity of the 1/f scale-invariant phenomenon has 
triggered in recent years the development of generic mecha- 
nisms describing complex systems, independent of their par- 
ticular context, in order to understand the "unifying" features 
of these systems [16-19]. 

To answer the question whether fluctuations in signals gen- 
erated by integrated physiological systems exhibit the same 
level of complexity, we analyze and compare the time series 
generated by two physiologic control systems under multiple- 
component integrated neural control — the human gait and 
the human heartbeat. We chose these two particular examples 
because human gait and heartbeat control share certain fun- 
damental properties, e.g., both originate in oscillatory centers. 



In the case of the heart, the pacemaker is located in the sinus 
node in the right atrium [20]. For gait, pacemakers called cen- 
tral pattern generators are thought to be located in the spinal 
cord [21]. 

However, these two systems are distinct, suggesting possi- 
ble dynamical differences in their output. For example, heart- 
beat fluctuations are primarily controlled by the involuntary 
(autonomic) nervous system. In contrast, while the sponta- 
neous walking rhythm is an automatic-like process, voluntary 
inputs play a major role. Further, gait control resides in the 
basal ganglia and related motor areas of the central nervous 
system, while the heartbeat is controlled by the sympathetic 
and parasympathetic branches of the autonomic nervous sys- 
tem [20, 22]. 

Previous studies show comparable two-point linear corre- 
lations and 1/f power spectra in heart rate [23-27] and hu- 
man gait [28-30], suggesting that differences in physiologic 
control may not be manifested in beat-to-beat and interstride 
interval fluctuations. Recent studies focusing on higher or- 
der correlations and nonlinear properties show that the human 
heartbeat exhibits not only 1/f fractal but also multifractal 
properties [31]. Since multifractal signals require many scal- 
ing indices to fully characterize their scaling properties, they 
may be considered to be more complex than those character- 
ized by a single fractal dimension, such as classical 1/f noise. 
Although the origins of the multifractal features in heartbeat 
dynamics are not yet understood, there is evidence that they 
relate to the complex intrinsic neuroautonomic regulation of 
the heart [31, 32]. Human gait, e.g., free unconstrained walk- 
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ing, is also a physiological process regulated by complex hi- 
erarchical feedback mechanisms involving supra-spinal inputs 
[21]. Moreover, recent findings indicate that the scaling prop- 
erties of gait fluctuations relate to neural centers on the higher 
supra-spinal level rather than to lower motor neurons or en- 
vironmental inputs [33, 34]. Thus it would be natural to hy- 
pothesize that the fluctuations in healthy unconstrained human 
gait exhibit similar fractal and multifractal features, and that 
human gait dynamics may belong to the same "complexity 
class" as cardiac dynamics. 

We employ two techniques — magnitude and sign decom- 
position analysis [35, 36] and multifractal analysis [37, 38] — 
to probe long-term nonlinear features, and to compare the lev- 
els of complexity in heartbeat and interstride interval fluctua- 
tions. To this end, we analyze interstride interval time series 
from 10 young healthy men (mean age 22 years) with no his- 
tory of neuromascular disorders [39]. Subjects walked contin- 
uously for 1 hour at a self-selected usual pace on level ground 
around a flat, obstacle-free, approximately oval, 400m long 
path. The interstride interval was measured using a ground re- 
action force sensor — ultra-thin force-sensitive switches were 
taped inside one shoe and data were recorded on an ambula- 
tory recorder using a previously validated method [40]. We 
compare the results of our gait analysis with results we have 
previously obtained [31, 35, 41, 42] from 6-hour long heart- 
beat interval records from 18 healthy individuals (13 female 
and 5 male, mean age 34 years) during daily activity (12:00 to 
18:00) [39]. 

As described below, we systematically compare the scal- 
ing properties of the fluctuations in human gait with those in 
the human heartbeat using power spectral analysis, detrended 
fluctuation analysis (DFA), magnitude and sign decomposi- 
tion analysis, and wavelet-based multifractal analysis, and we 
quantify linear and nonlinear features in the data over a range 
of time scales. 



II. METHODS 



A. Detrended fluctuation analysis (DFA) 



(if) We divide the profile y(k) into non-overlapping seg- 
ments of equal length n. 

(Hi) In each segment of length n, we fit y(k), using a poly- 
nomial function of order £ which represents the polynomial 
trend in that segment. The y coordinate of the fit line in each 
segment is denoted by y n (k). Since we use a polynomial fit 
of order £, we denote the algorithm as DFA-£. 

(z'v) The profile function y(k) is detrended by subtracting 
the local trend y n (k) in each segment of length n. In DFA- 
£, trends of order £ — 1 in the original signal are eliminated. 
Thus, comparison of the results for different orders of DFA-£ 
allows us to estimate the type of polynomial trends in the time 
series s(i). 

(v) For a given segment of length n, the root-mean-square 
(r.m.s.) fluctuation for this integrated and detrended signal 
s(i) is calculated: 



F(n) 
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(vz) Since we are interested in how F(n) depends on the 
segment length, the above computation is repeated for a broad 
range of scales n. 

A power-law relation between the average root-mean- 
square fluctuation function F(n) and the segment length n 
indicates the presence of scaling: 



F(n) 



(2) 



Thus, the DFA method can quantify the temporal organization 
of the fluctuations in a given signal s(i) by a single scaling 
exponent a — a self-similarity parameter which represents the 
long-range power-law correlation properties of the signal. If 
a = 0.5, there is no correlation and the signal is uncorrelated 
(white noise); if a < 0.5, the signal is anti-correlated; if a > 
0.5, the signal is correlated. The larger the value of a, the 
stronger the correlations in the signal. 

For stationary signals with scale-invariant temporal organi- 
zation, F(n) is related to the Fourier power spectrum S(f) 
and to the autocorrelation function C(n). For such signals, 



The DFA method was developed because conventional 
fluctuation analyses, such as power spectral, R/S and Hurst 
analysis cannot be reliably used to study nonstationary data 
[43, 45, 46]. One advantage of the DFA method is that it 
allows the detection of long-range power-law correlations in 
noisy signals with embedded polynomial trends that can mask 
the true correlations in the fluctuations of a signal. The DFA 
method has been successfully applied to a wide range of re- 
search fields in physics [47-51], biology [43, 52-55], and 
physiology [56-59]. 

The DFA method involves the following steps [44]: 
(z) Given the original signal s(i), where i = 1, ..,N max 
and N max is the length of the signal, we first form the profile 

function y(k) = 2~^i=i[ s (*) ~ ( s )]> where (s) is the mean. 
One can consider the profile y(k) as the position of a random 
walk in one dimension after k steps. 



S(f)~r p , where [ f3 = 2a - 1 ] (3) 

and a is the DFA scaling exponent (Eq. 2) [43, 44]. Thus sig- 
nals with 1/ / scaling in the power spectrum (i.e. f3 — 1) are 
characterized by DFA exponent a = 1. If 0.5 < a < 1, the 
correlation exponent 7 describes the decay of the autocorrela- 
tion function [43] : 

C(n) = (s(i)s(i + n)) ~ n~ 7 , where [ 7 = 2 - 2a }. (4) 

B. Magnitude and sign decomposition method 

Fluctuations in the dynamical output of physical and physi- 
ological systems can be characterized by their magnitude (ab- 
solute value) and their direction (sign). These two quantities 
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reflect the underlying interactions in a given system — the 
resulting "force" of these interactions at each moment deter- 
mines the magnitude and the direction of the fluctuations. Re- 
cent studies have shown that signals with identical long-range 
correlations can differ in the time organization of the magni- 
tude and sign of the fluctuations [35]. To assess the informa- 
tion contained in these fluctuations, the magnitude and sign 
decomposition method was introduced [35, 36]. This method 
involves the following steps: 

(z) Given the original signal s(i) we generate the increment 
series, As(i) = s(i + 1) — s(i). 

(ii) We decompose the increment series into a magnitude 
series |As(i)| and a sign series sgn(As(i)). 

(Hi) To avoid artificial trends we subtract from the magni- 
tude and sign series their average. 

(zv) We then integrate both magnitude and sign series, be- 
cause of limitations in the accuracy of the detrended fluctua- 
tion analysis method (DFA) for estimating the scaling expo- 
nents of anticorrelated signals (a < 0.5). 

(v) We perform a scaling analysis using 2nd order detrended 
fluctuation analysis (DFA-2) on the integrated magnitude and 
sign series. 

(vz) To obtain the scaling exponents for the magnitude and 
sign series we measure the slope of F(n)/n on a log-log plot, 
where F(n) is the root-mean-square fluctuation function ob- 
tained using DFA-2, and n is the scale. 

Fluctuations following an identical 1// scaling law can ex- 
hibit different types of correlations for the magnitude and the 
sign — e.g., a signal with anticorrelated fluctuations can ex- 
hibit positive correlations in the magnitude. Positive correla- 
tions in the magnitude series indicate that an increment with 
large magnitude is more likely to be followed by an increment 
with large magnitude. Anticorrelations in the sign series in- 
dicate that a positive increment in the original signal is more 
likely to be followed by a negative increment. Further, pos- 
itive power-law correlations in the magnitude series indicate 
the presence of long-term nonlinear features in the original 
signal, and relate to the width of multifractal spectrum [36] . 
In contrast the sign series relates to the linear properties of 
the original signal [36]. The magnitude and sign decomposi- 
tion method is suitable to probe nonlinear properties in short 
nonstationary signals, such as 1-hour interstride interval time 
series. 



C. Wavelet-based multifractal analysis 

Previously, analyses of the fractal properties of physio- 
logic fluctuations revealed that the behavior of healthy, free- 
running physiologic systems may often be characterized as 
l//-like [19, 23-27, 29, 34, 40, 60-71]. Monofractal signals 
(such as classical 1// noise) are homogeneous, i.e., they have 
the same scaling properties throughout the entire signal [72- 
76]. Monofractal signals can therefore be indexed by a single 
exponent: the Hurst exponent H [77]. 

On the other hand, multifractal signals are nonlinear and in- 
homogeneous with local properties changing with time. Mul- 
tifractal signals can be decomposed into many subsets charac- 



terized by different local Hurst exponents h, which quantify 
the local singular behavior and relate to the local scaling of the 
time series. Thus, multifractal signals require many exponents 
to fully characterize their properties [78-80]. The multifractal 
approach, a concept introduced in the context of multi-affine 
functions [81-84], has the potential to describe a wide class 
of signals more complex than those characterized by a single 
fractal dimension. 

The singular behavior of a signal s(t) at time to — 
\s(t) - P n {t)\ ~ \t - t |' l(to) for t -» t — is characterized 
by the local Hurst exponent h(to) where n < h(to) < n + 1 
and P n {t) is a polynomial fit of order n. To avoid an ad 
hoc choice of the range of time scales over which the local 
Hurst exponent h is estimated, and to filter out possible poly- 
nomial trends in the data which can mask local singularities, 
we implement a wavelet-based algorithm [38]. Wavelets are 
designed to probe time series over a broad range of scales 
and have recently been successfully used in the analysis of 
physiological signals [85-93]. In particular, recent stud- 
ies have shown that the wavelet decomposition reveals a ro- 
bust self-similar hierarchical organization in heartbeat fluc- 
tuations, with bifurcations propagating from large to small 
scales [42, 94, 95]. To quantify hierarchical cascades in gait 
dynamics and to avoid inherent numerical instability in the es- 
timate of the local Hurst exponent, we employ a "mean-field" 
approach — a concept introduced in statistical physics [1] — 
which allows us to probe the collective behavior of local sin- 
gularities throughout an entire signal and over a broad range 
of time scales. 

We study the multifractal properties of interstride interval 
time series by applying the wavelet transform modulus max- 
ima (WTMM) method [37, 38, 96] that has been proposed as a 
mean-field generalized multifractal formalism for fractal sig- 
nals. We first obtain the wavelet coefficient at time to from 
the continuous wavelet transform defined as: 

N 

W a (t ) = a- 1 ^s(t)^((t-t )/a), (5) 
t-i 

where s(t) is the analyzed time series, ip is the analyzing 
wavelet function, a is the wavelet scale (i.e., time scale of 
the analysis), and N is the number of data points in the time 
series. For tp we use the third derivative of the Gaussian, thus 
filtering out up to second order polynomial trends in the data. 
We then choose the modulus of the wavelet coefficients at 
each point t in the time series for a fixed wavelet scale a. 
Next, we estimate the partition function 

Z q (a) = Y^\W a (t)\\ (6) 

i 

where the sum is only over the maxima values of | W a (t ) | , and 
the powers q take on real values. By not summing over the en- 
tire set of wavelet transform coefficients along the time series 
at a given scale a but only over the wavelet transform modu- 
lus maxima, we focus on the fractal structure of the temporal 
organization of the singularities in the signal [96]. 



4 



We repeat the procedure for different values of the wavelet 
scale a to estimate the scaling behavior 



Z q (a) 



r(9) 



(7) 



In analogy with what occurs in scale-free physical systems, 
in which phenomena controlled by the same mechanism over 
multiple time scales are characterized by scale-independent 
measures, we assume that the scale-independent measures, 
r(q), depend only on the underlying mechanism control- 
ling the system. Thus by studying the scaling behavior of 
Z(a, q) ~ a T ^ we may obtain information about the self- 
similar (fractal) properties of the mechanism underlying gait 
control. 

For certain values of the powers q, the exponents r(q) have 
familiar meanings. In particular, r(2) is related to the scal- 
ing exponent of the Fourier power spectra, S(f) ~ I//' 3 , as 
(3 = 2 + t(2) [38]. For positive q, Z q (a) reflects the scaling 
of the large fluctuations and strong singularities in the signal, 
while for negative q, Z q (a) reflects the scaling of the small 
fluctuations and weak singularities [73, 80, 97]. Thus, the 
scaling exponents r(q) can reveal different aspects of the un- 
derlying dynamics. 

In the framework of this wavelet-based multifractal formal- 
ism, r(q) is the Legendre transform of the singularity spec- 
trum D(h) defined as the Hausdorff dimension of the set of 
points t in the signal s(t) where the local Hurst exponent is 
h. Homogeneous monofractal signals — i.e., signals with a 
single local Hurst exponent h — are characterized by linear 
r(q) spectrum: 



r(q) =qH-l, 



(8) 



where H = h = dr(q)/dq is the global Hurst exponent. On 
the contrary, a nonlinear r(q) curve is the signature of nonho- 
mogeneous signals that display multifractal properties — i.e., 
h(t) is a varying quantity that depends upon t. 



III. RESULTS 

In Fig. 1 we show two example time series: (i) an inter- 
stride interval time series from a typical healthy subject dur- 
ing w 1 hour (3, 000 steps) of unconstrained normal walking 
on a level, obstacle-free surface (Fig. la) [39]; (ii) consecu- 
tive heartbeat intervals from w 1 hour (3, 000 beats) record of 
a typical healthy subject during daily activity (Fig. lb) [39]. 
Both time series exhibit irregular fluctuations and nonstation- 
ary behavior characterized by different local trends; in fact it 
is difficult to differentiate between the two time series by vi- 
sual inspection. 

We first examine the two-point correlations and scale- 
invariant behavior of the time series shown in Fig. 1. Power 
spectra S(f) of the gait and heartbeat time series (Fig. 2a) 
indicate that both processes are described by a power-law re- 
lation S(f) ~ 1/ f 13 over more than 2 decades, with exponent 
(3 w 1. This scaling behavior indicates self-similar (fractal) 
properties of the data suggestive of an identical level of com- 
plexity as quantified by this linear measure. We obtain similar 



results for the interstride interval times series from all subjects 
in our gait database: j3 = 0.9 ± 0.08 (group mean ± std. dev.) 
in agreement with previous results [34]. 



A. Detrended fluctuation analysis (DFA) 

Next, to quantify the degree of correlation in the interstride 
and heartbeat fluctuations we apply the DFA method, which 
also provides a linear measure: plots of the root-mean-square 
fluctuation function F(n) vs. time scale n (measured in stride 
or beat number) from a second-order DFA analysis (DFA-2) 
[44—46] indicate the presence of long-range power-law corre- 
lations in both gait and heartbeat fluctuations (Fig. 2b). The 
scaling exponent a w 0.95 for the heartbeat signal is very 
close to the exponent a w 0.9 for the interstride interval sig- 
nal, estimated over the scaling range 6 < n < 600. We ob- 
tain similar results for the remaining subjects: a = 0.87±0.03 
(group mean ± std. dev.) for the gait data and a = 1.04±0.08 
for the heartbeat data, in agreement with [34]. 

The results of both power spectral analysis and the DFA 
method indicate that gait and heartbeat time series have sim- 
ilar scale-invariant properties suggesting parallels in the un- 
derlying mechanisms of neural regulation. 



B. Magnitude and sign decomposition method 

To probe for long-term nonlinear features in the dynam- 
ics of interstride intervals we employ the magnitude and 
sign decomposition analysis [35, 36]. Previous studies have 
demonstrated that information about the nonlinear properties 
of heartbeat dynamics can be quantified by long-range power- 
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FIG. 1: Representative records of (a) interstride interval time series 
from a healthy subject and (b) consecutive heartbeat intervals from a 
healthy subject 
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FIG. 2: (a) Power spectra of the gait series (A) and heartbeat series 
(•) displayed in Fig. 1. Plots of the root-mean-square fluctuation 
function F(n) vs. time scale n (measured in stride or beat number) 
from second-order DFA-2 analysis for (b) the interstride and heart- 
beat interval time series indicating similar power-law correlations, 
and (c) the magnitude series of the interstride and heartbeat incre- 
ments showing a surprising difference in the nonlinear properties of 
the two time series. 



law correlations in the magnitude of the increments in heart- 
beat intervals [35]. Further, correlations in the magnitude are 
associated with nonlinear features in the underlying dynam- 
ics, while linear signals are characterized by an absence of 
correlations (random behavior) in the magnitude series. To 
quantify the correlations in the magnitude of the interstride 
increments we apply the DFA-2 method to the data displayed 



in Fig. la. Our results show that the magnitude series of the 
interstride increments exhibits close to random behavior with 
correlation exponent a « 0.5 (denoted by (A) in Fig. 2c). 
In contrast, the magnitude series of the heartbeat increments 
(Fig. lb) exhibits strong positive correlations over more than 
two decades characterized by exponent a = 0.7 (denoted 
by (•) in Fig. 2c). A surrogate test [98, 99] eliminating the 
nonlinearity in the heartbeat time series by randomizing the 
Fourier phases but preserving the power spectrum leads to ran- 
dom behavior {a = 0.5) in the magnitude series [35]. Thus 
the striking difference in the magnitude correlations of gait 
and heartbeat dynamics (both of which are under multilevel 
neural control) raises the possibility that these two physiologic 
processes belong to different classes of complexity whereby 
the neural regulation of the heartbeat is inherently more non- 
linear, over a range of time scales, than the neural mechanism 
of gait control. Our observation of a low degree of nonlin- 
earity in the gait time series is supported by the remaining 
subjects in the group: over time scales 6 < n < 600, we 
obtain exponent a mag = 0.57±0.03 (group mean ± std. dev.) 
for the gait time series, which is significantly lower than the 
corresponding exponent a mag = 0.75 ± 0.06 obtained for the 
heartbeat data (p = 4.8 x 10~ 6 , by the Student's t-test). 



C. Wavelet-based multifractal analysis 

To further test the long-term nonlinear features in gait dy- 
namics we study the multifractal properties of interstride time 
series. We apply the Wavelet Transform Modulus Maxima 
(WTMM) method [38, 96] — a "mean-field" type approach 
to quantify the fractal organization of singularities in the sig- 
nal. We characterize the multifractal properties of a signal 
over a broad range of time scales by the multifractal spectrum 

We first examine the time series shown in Fig. 1. For the 
gait time series, we obtain a r(q) spectrum which is practi- 
cally a linear function of the moment q, suggesting that the 
gait dynamics exhibit monofractal properties (Fig. 3a). This 
is in contrast with the nonlinear r(q) spectrum for the heart- 
beat signal (Fig. 3a) which is indicative of multifractal behav- 
ior [37, 38]. Further when analyzing the remaining interstride 
interval recordings we find close to linear r(q) spectra for all 
subjects in the gait group (Fig. 3b). Calculating the group av- 
eraged r(q) spectra we find clear differences: multifractal be- 
havior for the heartbeat dynamics and practically monofractal 
behavior for the gait dynamics (Fig. 3c). Specifically we find 
significant differences between the gait and heartbeat r(q) 
spectra for negative values of the moment q\ for positive val- 
ues of q, the scaling exponents r(q) take on similar values. 
This is in agreement with the similarity in power spectral and 
DFA scaling exponents for gait and heartbeat data, which cor- 
respond to r(q = 2) (Fig. 2). However, the heartbeart r(q) 
spectrum is visibly more curved for all moments q compared 
with the gait r(q) spectrum which may be approximately fit 
by a straight line, indicative of a low degree of nonlinearity 
in the interstride time series. Thus our results show consistent 
differences between the nonlinear and multifractal properties 
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FIG. 3: Multifractal analysis: (a) Multifractal spectrum r(q) for the 
individual records shown in Fig. 1, where r is a scaling index asso- 
ciated with different moments q (Eq. 7). A monofractal signal corre- 
sponds to a straight line for r(q), while for multifractal signals r(q) 
is a nonlinear function of q. The values of r(q = 2) for both gait 
and heartbeat time series are very close, in agreement with our find- 
ings based on DFA-2 correlation analysis (Fig. 2b). (b) Multifractral 
spectra r(q) for all ten subjects in our database [39] exhibit close to 
linear dependence on the moment q, suggesting monofractal behav- 
ior, in contrast to the nonlinear r(q) spectra reported for heartbeat 
recordings [102]. (c) Group average multifractal spectra r(q) for the 
gait and heartbeat subjects in our database [39]. The results show a 
consistent monofractal (almost linear) behavior for the gait time se- 
ries, in contrast with the multifractal behavior of the heartbeat data. 



of gait and heartbeat time series. 

Previous studies have shown that reducing the level of phys- 
ical activity under a constant routine protocol does not change 
the multifractal features of heartbeat dynamics, while block- 
ing the sympathetic or parasympathetic tone of the neuro- 
autonomic regulation of the heart dramatically changes the 
multifractal spectrum, thus suggesting that the observed fea- 
tures in cardiac dynamics arise from the intrinsic mechanisms 
of control [32]. Similarly, by eliminating polynomial trends in 
the interstride interval time series corresponding to changes 
in the gait pace using DFA and wavelet analyses, we find 
scaling features which remain invariant among individuals. 
Therefore, since different individuals experience different ex- 
trinsic factors, the observed lower degree of nonlinearity as 
measured by the magnitude scaling exponent and the close-to- 
monofractal behaviour characterized by practically linear r(q) 
spectrum appear to be a result of the intrinsic mechanisms of 
gait regulation. These observations suggest that while both 
gait and heartbeat dynamics arise from layers of neural con- 
trol with multiple component interactions, and exhibit tem- 
poral organization over multiple time scales, they nonetheless 
belong to different complexity classes. While both gait and 
heartbeat dynamics may be a result of competing inputs inter- 
acting through multiple feedback loops, differences in the na- 
ture of these interactions may be imprinted in their nonlinear 
and multifractal features: our findings suggest that while these 
interactions in heartbeat dynamics are of a nonlinear character 
and are represented by Fourier phase interactions encoded in 
the magnitude scaling and the multifractal spectrum, feedback 
mechanisms of gait dynamics lead to decreased interactions 
among the Fourier phases. 



D. Further validation of gait results 

These new findings are supported by our analysis of a sec- 
ond group of gait subjects. We analyze interstride intervals 
from an additional group of 7 young healthy subjects (6 male, 
1 female, mean age 28 years) recorded using a portable ac- 
celerometer [100]. Subjects walked continuously for w 1 hour 
at a self-selected pace on an unconstrained outdoor walking 
track in a park environment allowing for slight changes in el- 
evation and obstacles related to pedestrian traffic. The stride 
interval time series in this case were obtained from peak-to- 
peak intervals in the accelerometer signal output in the direc- 
tion of the subjects' vertical axis [101]. Compatibility of the 
ground reaction force sensor used for the gait recordings of 
the first group [40] with the accelerometer device, and strong 
correlation between outputs of the two devices was reported 
inRef. [100]. 

We find that for this second group the two-point correlation 
exponent a, as measured by the DFA method a = 0.90 ±0.1 
(group mean ± std. dev.) is similar to the group aver- 
age exponent of the first gait group (a = 0.87 ± 0.03) and 
also the heartbeat data (a = 1.04 ± 0.08). In contrast, 
we find again a significantly lower degree of nonlinearity, as 



measured by the magnitude exponent a„ 



0.62 ± 0.04 



and the r{q) spectrum, compared with heartbeat dynamics 
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a mag = 0.75 ± 0.06 (p = 1.3 x 10~ 3 , by the Student's t- 
test) (Fig. 2c and Fig. 3c). On the other hand, the group aver- 
aged value of a mag is slightly higher compared with the first 
gait group (a mag = 0.57 ± 0.03), and this is associated with 
slightly stronger curvature in the r(q) spectrum for the second 
gait group. This may be attributed to the fact that the second 
group walked in a natural park environment where obstacles, 
changes in elevation and pedestrian traffic may possibly re- 
quire the activation of higher neural centers. 

The present results are related to a physiologically-based 
model of gait control where specific interactions between neu- 
ral centers are considered [12, 13]. In this model a lower 
degree of nonlinearity (and close-to-linear monofractal r(q) 
spectrum) reflects increased connectivity between neural cen- 
ters, typically associated with maturation of gait dynamics in 
adults. The present results are also consistent with studies that 
used a different approach to quantify the dynamics of giat, 
based on estimates of the local Hurst exponents, and reported 
only weak multifractality in gait dynamics [14, 15]. 

IV. SUMMARY 

In summary, we find that while the fluctuations in the out- 
put of both gait and heartbeat processes are characterized by 



similar two-point correlation properties and l//-like spectra, 
they belong to different classes of complexity — human gait 
fluctuations exhibit linear and close to monofractal properties 
characterized by a single scaling exponent, while heartbeat 
fluctuations exhibit nonlinear multifractal properties which in 
physical systems have been connected with turbulence and re- 
lated multiscale phenomena [37, 83, 84, 103]. 

These findings are of interest because they underscore the 
limitations of traditional two-point correlation methods in 
characterizing physiologic and physical time series. In ad- 
dition, these results suggest that feedback on multiple time 
scales is not sufficient to explain different types of 1/ f scal- 
ing and scale-invariance, and highlight the need for the devel- 
opment of new models [104-107] that could account for the 
scale-invariant outputs of different types of feedback systems. 
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